function A = subtable(f,fexact)
precision = 0;
gi=1;
while(precision <7)
    g = gauss(f,gi);
    grel = (g-fexact)/fexact;
    precision = -log10(grel);
    gi=gi+1;
end
precision = 0;
ci=1;
while(precision <7)
    c = clenshaw_curtis(f,ci);
    crel = (c-fexact)/fexact;
    precision = -log10(crel);
    ci=ci+1;
end
precision = 0;
ri=1;
reval=0;
while(precision <7)
    r = romberg(f,ri);
    rrel = (r-fexact)/fexact;
    precision = -log10(rrel)
    ri=ri*2;
end
A=[gi ci ri];